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^v/l ' Abstract. We consider Jacobi matrices with zero diagonal and off-diagonals 

' given by elements of the hull of the Fibonacci sequence and show that the 

spectrum has zero Lebesgue measure and all spectral measures are purely 
singular continuous. In addition, if the two hopping parameters are distinct 
but sufficiently close to each other, we show that the spectrum is a dynamically 
defined Cantor set, which has a variety of consequences for its local and global 
fractal dimension. 
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1. Introduction 



C^ , The study of the Fibonacci operator was initiated in the early 1980's by Kohmoto 

et al. [17] and Ostlund et al. [23|. At that point in time, the interest in this model 
was based mainly on the existence of an exact renormalization group procedure 
and the appearance of critical eigenstates and zero measure Cantor spectrum. Only 
shortly thereafter, Shechtman et al. [28' reported their discovery of structures, now 
'^b ' called quasicrystals, whose diffraction exhibits sharp Bragg peaks with rotational 

^\J , symmetries that ensure the aperiodicity of the structure in question. 

^D ' The Fibonacci sequence is the central model of a quasicrystal in one dimen- 

^^ , sion. Indeed, it is aperiodic and pure point diffractive and, moreover, it belongs to 

t^^ ' virtually all classes of mathematical models of quasicrystals that have since been 

C3 \ proposed. We refer to the reader to T for a relatively recent account of the math- 

ematics related to the modelling and study of quasicrystals. 

Thus, the study of the Fibonacci operator was further motivated by the inter- 
est in electronic spectra and transport properties of one-dimensional quasicrystals. 
Jv>( \ Consequently, apart from the almost Mathieu operator, the Fibonacci operator has 

5^ ' been the most heavily studied quasi-periodic operator in the last three decades; 

5^ \ compare the survey articles |51 [3D] . 

Partly due to the choice of the model in the foundational papers pTl [23] , the 
mathematical literature on the Fibonacci operator has so far only considered the 
diagonal model, that is, a discrete one-dimensional Schrodinger operator with po- 
tential given by the Fibonacci sequence. Given the connection to quasicrystals 
and hence aperiodic point sets, and particularly cut-and-project sets, it is however 
equally (if not more) natural to study the off-diagonal model. That is, one con- 
siders a Jacobi matrix with zero diagonal and off-diagonals given by the Fibonacci 
sequence. This would correspond to hopping over unequal distances, the sites being 
points on the real line generated by the cut-and-project scheme that produces the 
Fibonacci point set. While we will not explain these concepts here in detail, we do 
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refer the reader to [1] for background and we will describe the resulting model care- 
fully in the next section. For further motivation to study the off-diagonal model, 
we mention that it has been the object of interest in a number of physics papers; 
see, for example, [HI [H [HI [lH [SH [32] . 

In fact, our interest in the off-diagonal model was triggered by the numeri- 
cal results for separable higher-dimensional models obtained recently by Even-Dar 
Mandel and Lifshitz in jTSl [M] . We intend to explore the observations reported in 
these papers and prove some of the phenomena rigorously, and for this task, we will 
need the basic spectral results for the off-diagonal Fibonacci operator obtained in 
this paper. 

The paper is organized as follows. In Section [2| we describe in detail the model 
we study and the results we obtain for it. The results concerning the spectrum as a 
set are proved in Section [3| and the results concerning the spectral type are proved 
in Section [D 

Acknowledgment. We wish to thank Serge Cantat for useful conversations. 



2. Model and Results 

In this section we describe the off-diagonal Fibonacci operator and the results 
for it that are proved in subsequent sections. Let a, b be two positive real numbers 
and consider the Fibonacci substitution, 

S{a) = ab, S{b) = a. 

This substitution rule extends to finite and one-sided infinite words by concate- 
nation. For example, S{aba) = abaab. Since S{a) begins with a, one obtains a 
one-sided infinite sequence that is invariant under S by iterating the substitution 
rule on a and taking a limit. Indeed, we have 

(1) S'^ia) = S''-\S{a)) = S''-^{ab) = S''-\a)S''-^ib) = S''-\a)S''-^{a). 

In particular, S^'^a) starts with iS''^~^(a) and hence there is a unique one-sided 
infinite sequence u, the so-called Fibonacci substitution sequence, that starts with 
S'^{a) for every k. The hull fla^b is then obtained by considering all two-sided 
infinite sequences that locally look like u, 

^a.b = {w G {a, b}^ : every subword of a; is a subword of u}. 

It is known that, conversely, every subword of u is a subword of every lu e ^a.b- I^^ 
this sense, u and all elements of the hull uj look exactly the same locally. 

We wish to single out a special element of fta.b- Notice that ba occurs in u and 
that S'^{a) = aba begins with a and 5^(6) = ab ends with b. Thus, iterating S^ on 
b\a, where | denotes the eventual origin, we obtain as a limit a two-sided infinite 
sequence which belongs to ^a,b and coincides with u to the right of the origin. This 
element of fla.b will be denoted by cus- 

Each uj G ila^b generates a Jacobi matrix H^ acting in ^^(Z), 

iHu;1p)n = UJn+l-ipn+l + UJ„1pn-l- 
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With respect to the standard orthonormal basis {Sn}nez of ^^(Z), where (5„ is one 
at n and vanishes otherwise, this operator has the matrix 
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and it is clearly self-adjoint. 

This family of operators, {Huj}ujen^ i,, is called the off-diagonal Fibonacci model. 
Of course, the structure of the Fibonacci sequence disappears when a — b. In this 
case, the hull consists of a single element, the constant two-sided infinite sequence 
taking the value a = b, and the spectrum and the spectral measures of the associated 
operator H^ are well understood. For this reason, we will below always assume that 
a ^ b. Nevertheless, the limiting case, where we fix a, say, and let b tend to a is of 
definite interest. 

Our first result concerns the spectrum of iJ^. For S' C M, we denote by dim/f S* 
the Hausdorff dimension of S and by dimes' the box counting dimension of S (which 
is then implicitly claimed to exist); see [8] or [24] for the definition of these fractal 
dimensions. 

Theorem 1. Suppose a,b > and a ^ b. Then, there exists a compact set T,a.b C M 
such that a{H^) = Y^a.b for every lu G ^a,b, arid 
(i) Sab has zero Lebesgue measure. 

(ii) The Hausdorff dimension diuiH Sa,fc is an analytic function of a and b. 

(iii) < dimH^a,b < 1- 

2 I .2 

Also, there exists Eq > such that if °„ , < 1 + Eq {in other words, if a and b 
are close enough), then 

(iv) The spectrum Tia.b is a Cantor set that depends continuously on a and b in 

the Hausdorff metric. 
(v) For every small 5 > Q and every E G Saf,, we have 

dimjf {{E -d,E + d)n E,^b) = dims {{E -d,E + d)n E^.b) 

= dims Sa.fc = dims Sa.b. 

(vi) Denote a — dim^ '^a,b, then the Hausdorff a-measure of T.a^b is positive 
and finite. 

Properties (ii)-(vi) are related to hyperbolicity of the trace map described in 
the next section, see [21 |4j |9] . It is reasonable to expect that in fact the properties 
(iv)-(vi) for Y.a.b hold for all positive a ^ b, but one of the steps in the proof is 
currently missing as will be explained in Section [31 

Next, we turn to the spectral type of H^i. 

Theorem 2. Suppose a, 6 > and a ^ b. Then, for every to G i^a.b, H^^ has purely 
singular continuous spectrum. 



T^(n,£;) = — I 2 n 
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Throughout the rest of the paper we will only consider a, 6 > with a ^ h. 
Theorem [1] is proved in Section [3] and Theorem [2] is proved in Section ID 

3. The Trace Map and its Relation to the Spectrum 

The spectral properties of H^ are closely related to the behavior of the solutions 
to the difference equation 

(2) UJn+lUn+l + LUnUn~l = EUn- 

Denote 

\UJnUn-l 

Then u solves ([2]) (for every rt G Z) if and only if U solves 

(3) [/„=r^(n,^)[/„_i, 
(for every n G Z), where 

}_fE -1' 

Note that deiT^{n,E) ^ 1. Iterating dH), we find 

[/„ = A'f^(n,^)f/o, 

where 

M^(n,^)=T^(n,^)x ■■■T^{\,E). 

With the Fibonacci numbers {-Ffc}, generated by Fq = F\ = 1, Ffe+i = i^^ + i^fe_i 
for fc > 1, we define 

Xk=Xk{E) = \^TM^XFk.E). 
For example, we have 

M..(F„i.)4(j -; 

\ (E -\\^(E -f\ 1 /S2_^2 _^ 

h\h^ Q)~a\o? )-^[ EP -h\ 

1 (E -l\ 1 /£; -1\ 1 (E -l\ _ 1 fE^'-Ea^-Eb^ -E^ + b^ 

a[a^ ) b[b^ ) a[a^ ) ^ a% [ E^a^ - a* -Ea 

and hence 

E E^-a^-b^ E^- 2Ea^ - Eb^ 

(') "^^2^' "^ = 2^6 ' "^^ 2^^6 ■ 

Lemma 1. We have 

(5) Xk+i = 2xkXk-i - Xk-2 
for k > 2. Moreover, the quantity 

(6) Ik = a^fe+i + xl + xl^i - 2xk+iXkXk-i - 1 
is independent of both k and E and it is given by 

. {a' + br , 
4a252 



Mu^XF^^E) - - ^ 2 nlTlL2 n/Tl„2 n/-T^l Z?2„2 „4 77„2 
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Proof. Since ujg restricted to {n > 1} coincides with u and the prefixes Sk of u of 
length Fk obey Sk+i — SkSk-i for fc > 2 by construction (cf. ^), the recursion ^ 
follows as in the diagonal case; compare [51 [SI [^. This recursion in turn implies 
readily that ([6]) is fc-independent. 

In particular, the x^-'s are again generated by the trace map 

T{x,y,z) = {2xy- z,x,y) 

and the preserved quantity is again 

I{x,y,z) = x +y +z —2xyz — l. 

The only difference between the diagonal and the off-diagonal model can be found 
in the initial conditions. How are xijXca;-! obtained? Observe that the trace 
map is invertible and hence we can apply its inverse twice to the already defined 
quantity {x3,X2,xi). We have 

T~^{x,y,z) = {y,z,2yz~x) 

and hence, using ([4]), 

(a;i,a;o,x-i) = T~'^{x3,X2,xi) 

_ 2 / E^-2Ea^~Eb^ E'^-a'- b^ E\ 

■ E^-a^-y" E (E^ - a^ - b^)E E^ ^ 2Ea^ ~ Eb^ 



T- 



2ab '2a' Aa^ 2a^ 

E^-a^-b^ E E' 



2ab ' 2a' 2b 

E_ E_ E^ E^-a^- b^ 
2a'26' Aab 2ab 

E E a^ + r 



2a' 26' 2a6 



It follows that 



I{xk+l,Xk,Xk^l) = I{xi,Xo,X^i) 



E^ , E^ , (a2+62)2 ^2(^2 



4a2 4fe2 4a252 8a2&2 

^ (a2 + b2)2 

4a262 ^ 

for every A: > 0. D 

It is of crucial importance for the spectral analysis that, as in the diagonal case, 
the invariant is energy-independent and strictly positive when a ^ b\ 

Lemma 2. The spectrum of H^ is independent of lu and may be denoted by T,a,b- 
With 

ak = {EeR: \xk\ < 1}, 
we have 

(7) Ea^b = (I CTfe U CTfe+i. 

k>l 
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Moreover, for every E e Tia,b and k >2, 

and for E 'Sa,b, \xk\ diverges super- exponentially. 

Proof. It is well known that the huU fla^ together with the standard shift transfor- 
mation is niininiaL In particular, every to S Qa,b may be approximated pointwise by 
a sequence of shifts of any other lu G ^a,b- The associated operators then converge 
strongly and we get a{H^) C a{Hcj). Reversing the roles of w and uj, the first claim 
follows. 

So let Y,a,b denote the common spectrum of the operators H^^, uj G ^a,b- We 
have \\H^\\ < max{2a,26}. Thus, Sa,6 C [- max{2a, 26}, max{2a, 26}] =: Ia,b- For 
E G Ia,b, we have that at least one of xi, xq belongs to [—1, 1]. This observation 
allows us to use the exact same arguments Siito used to prove ([7]) for the diagonal 
model in ^5] , 

The only point where care needs to be taken is the claim that (Jk is the spectrum 
of the periodic Jacobi matrix obtained by repeating the values lUs takes on {1 < n < 
Fk} periodically on the off-diagonals. This, however, follows from the general theory 
of periodic Jacobi matrices, which relies on the diagonalization of the monodromy 
matrix (which is Afc^^ [Fk , E) in this case) in order to obtain Floquet solutions and 
in particular discriminate between those energies that permit exponentially growing 
solutions and those that do not. This distinction works just as well here, but one 
needs to use that the w„ 's that enter in the C/„ 's are uniformly bounded away from 
zero and infinity. 

Thus, after paying attention to this fact, we may now proceed along the lines of 
Siito. Let us describe the main steps of the argument. Since at least one of Xi, Xq 
belongs to [—1,1], we have a result analogous to [5S1 Lemma 2] with the same proof 
as given there. Namely, the sequence {xk}k>Q is unbounded if and only if there 
exists k such that |xfc| > 1 and |a::/c+i| > 1. Moreover, we then have \xk+i\ > c^' for 
some c > 1 and alH > 0. This shows 

CTfe Uo-fc+i = [J (Jk+l- 
l>0 

Using now the fact that the F^. periodic Jacobi matrices with spectrum ak converge 
strongly to H^^^ , we obtain 



k>ll>0 k>l k>l 

since the spectra ak and (Jk+i are closed sets. Thus, we have one inclusion in ([7]). 
Next, suppose E G nfe>i "■*: Ucfe+i If fc > 1 is such that \xk\ > 1, then l^fc^ij < 1 
and |xfe4-i| < 1. Since we have 



4-^1 +^1+ 4-1 - 2xk+iXkXk-i - 1 = ^^2j,2 1' 



this implies 



xk = xk+ixk-i ± ( 1 - 4+1 - 4-1 + 4+i4-i 



,2 I l2\2 \ 1/2 



4a262 
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and hence 

\xk\ < \xk+iXk-i\ + ( (1 - 4+i)(l - -^k-l) 



4a262 



which, using |a;fc_i| < 1 and |a:fe+i| < 1 again, imphes the estimate ([8]) for E £ 
nfe>i ^fc U o'fe+i- We will show in the next section that the boundedness of the 
sequence {xk}k>o implies that, for arbitrary to £ Oa.fc, no solution of the difference 
equation ^ is square-summable at -l-cxo. Consequently, such E^s belong to Sa.fcO 
This shows the other inclusion in ([7]) and hence establishes it. Moreover, it follows 
that ^ holds for every E e Y.a,b- 

Finally, from the representation ([7]) of Y.a.b and our observation above about 
unbounded sequences {xk}k>o, we find that \xk\ diverges super-exponentially for 
E ^ Sa,fc- This concludes the proof of the lemma. D 

Lemma 3. For every E ^M., there is ^{E) > such that 

hm -log\\MUn,E)\\^j{E), 
uniformly in lu € ^a.b- 
Proof. This follows directly from the uniform subadditive ergodic theorem; compare 

[nmansiEni. □ 

Lemma 4. The set Z^^h '■— {E e M : ^{E) = 0} has zero Lebesgue measure. 

Proof. This is one of the central results of Kotani theory; see fT9| and also [7|. 
Note that these papers only discuss the diagonal model. Kotani theory for Jacobi 
matrices is discussed in Carmona-Lacroix j2j and the result needed can be deduced 
from what is presented there. For a recent reference that states a result sufficient 
for our purpose explicitly, see Remling [26] . D 

Lemma 5. We have I]a,b = Za.b. 

Proof. The inclusion T,a.b 3 Za.b holds by general principles. For example, one can 
construct Weyl sequences by truncation when ^{E) = 0. The inclusion S^.h ^ Za.b 
can be proved in two ways. Either one uses the boundedness of Xk for energies 
E G Sa,fc to prove explicit polynomial upper bounds for ||A/i^(n, _E)|| (as in [16j for 
uj = ujs ov in TT] for general uj € ^a,b), or one combines the proof of the absence of 
decaying solutions at +oo for E G Sat given in the next section with Osceledec's 
Theorem, which states that "f{E) > would imply the existence of an exponentially 
decaying solution at +oo. Here we use one more time that C/„ is comparable in 
norm to (u„,u„„i)^. D 

Proof of Theorem\^ The existence of the uniform spectrum Y,a,b was shown in 
Lemma[3]and the fact that E^.h has zero Lebesgue measure follows from Lemmas|3] 
and El The set of bounded orbits of the restriction of the trace map T : R'^ ^ M"^ to 

the invariant surface /(x, y,z) = C = ^ 4^2 ^2 1, C > 0, is hyperbolic; see [3J (and 

also [9] for C sufhciently small and j4j for C sufficiently large). Due to Lemma [21 



This follows by a standard argument: If E ^ ^a,bt then (H^ — E)~^ exists and hence 
{H^ — E)^^8o is an £^(Z) vector that solves ^ away from the origin. Choosing its values for 
n > 1, say, and then using ^ to extend it to all of Z, we obtain a solution that is square-summable 

at +00. 
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the points of the spectrum correspond to the points of the intersection of the hne 
of the initial conditions 

, _(/E E a^ + b^\ ^ 

''''''= U2^' 26' ^^rj^^^ 

with the stable manifolds of the hyperbolic set of bounded orbits. Properties (ii) 
and (iii) can be proved in exactly the same way as Theorem 6.5 in [3]. The line la,b 
intersects the stable lamination of the hyperbolic set transversally for sufficiently 
small C > 0, as can be shown in the same way as for the diagonal Fibonacci 
Hamiltonian with a small coupling constant; see [9]. Therefore the spectrum J^a,b 
for close enough a and 6 is a dynamically defined Cantor set, and the properties 
(iv)-(vi) follow; see [H [HI [U [H [HI [21] and references therein. D 

Notice that a proof of the transversality of the line Za.fc to the stable lamination of 
the hyperbolic set of bounded orbits for arbitrary a ^ b would imply the properties 
(iv)-(vi) for these values of a and b. 

4. Singular Continuous Spectrum 

In this section we prove Theorem[2l Given the results from the previous section, 
we can follow the proofs from the diagonal case quite closely. 

Proof of Theorem [^ Since the absence of absolutely continuous spectrum follows 
from zero measure spectrum, we only need to show the absence of point spectrum. 
It was shown by Damanik and Lenz lOj that, given any lu £ ^a,b and fc > 1, the 
restriction oi lu to {n > 1} begins with a square 

Wl . . . U!2Fk ■ ■ ■ — UJl . . . UJp^UJl . . . UJPk ■ ■ ■ 

such that (jji . . . wfj. is a cyclic permutation of S^{a). By cyclic invariance of the 
trace, it follows that TyM^{Fu,E) = 2xk{E) for every E. 

The Cayley- Hamilton Theorem, applied to AI^{Fk,E), says that 

M^{Fk,Ef - (TrM^(Ffe, S)) M^{Fk,E) + 1 = 0, 

which, by the observations above, translates to 

M^{2Fk,E) - 2xkMUFk,E) + 1 = 0. 

Ii E £ 'Ea.b and u is a solution of the difference equation ^, it therefore follows 
that 

U{2Fk + 1) - 2xkU{Fk + 1) + C/(l) = 0. 

If u does not vanish identically, this shows that u„ y^- as n — > oo since the Xk 's are 
bounded above and the a;„'s are bounded below away from zero. In particular, if 
E G Sa,fc, then no non-trivial solution of ([2]) is square-summable at +cxd and hence 
E is not an eigenvalue. It follows that the point spectrum of H^ is empty. D 
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